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<N ■ 1. INTRODUCTION 

O 

Q . The most general formulation of phenomenological models of magnets (or spin sys- 

tems) which includes all the known completely integrable ones, has the following form: 

00 : (1-1) 

\ + Ot^Uyy = -Ro(S, S^, Sy), 

, where S = S{x,y,t) is the magnetization vector, Fo(,) is a vector-function, u = u{x,y) 

is an auxiliary field, Rq{, ) is a scalar function, J is a set of constants characterizing the 
magnet, and = ±1. 
\0 ' The function Fq usually takes the form: 



> 
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o 



SA^^ + Fi, (1.2) 



where F^jf is the functional of the crystal's free energy (throughout the paper the symbol 
6/6 stands for the variational derivative). The first term in the right hand side was 
suggested by Landau and Lifshits [1] to describe the exchange interactions. 
. The representation (1. !)-(!. 2) is often inconvenient for solving problems. One would 

! like to deal with more tractable forms of the equations (1.1), which, in turn, requires 

5^ \ introduction of new dependent variables. Apparently, such a variable, the stereographic 

projection, has been used for the first time in paper [2] to describe the instanton solutions 
in the two-dimensional 0(3) cr-model (the 2D stationary Heisenberg ferromagnet). Later 
it was exploited in various situations, see, e. g. [3-5]. 

In the present paper we show on examples of three models - the deformed Heisen- 
berg, the Landau - Lifshits, and the Ishimori magnets - that it is helpful to introduce 
the corresponding canonical variables. In particular, they allow to simplify significantly 
certain calculations, as compared to the usage of the S variable, and, more substantially, 
to clarify a set of questions important both from physical and mathematical viewpoints. 
Another argument in their favor is that in these variables the models fit in a class of 
models admitting a differential-geometric interpretation intensively studied recently [6]. 

The model of deformed Heisenberg magnet was suggested in [7] where also an exact 
solution of it for the case of trivial background was obtained by the inverse scattering 
method, and the conservation laws were calculated. In doing so, it was shown that 
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perturbations localized in the space are spreaded, that is, the solutions are instable. The 
gauge equivalence of this model and the nonlinear Schrodinger equation with an integral 
nonlinearity was established in [8] and [9]. The matrix Darboux transform method was 
applied in the paper [10], where exact solutions of the model where calculated on the 
background of new spiral-logarithmic structures. 

The Landau - Lifshits equation is a subject of vast studies. In particular, the Lax rep- 
resentation and conservation laws for it in the completely anisotropic case have first been 
obtained in [11], soliton solutions were found by the dressing method in [12]. Hamilton 
aspects of the equation were analyzed in detail in recent paper [13]. 

The Ishimori magnet was also considered in many papers. In particular, series of 
exact solutions were obtained in [14] by the inverse scattering and d - dressing methods, 
the Darboux transform was applied to it in [15] and [16] (in [16] - on the background of 
spiral structures). Notice also the important paper [5], where the gauge equivalence of 
the Ishimori-II and Davy-Stewartson-II models was established. 

The structure of this paper is as follows. In section 2 we consider the deformed Heisen- 
berg magnet model, define the canonical variables and analyze stability of the solutions. 
In section 3 the Landau - Lifshits equation is obtained in terms of the stereographic 
projection and a stationary version of this equation is studied. Finally, in section 4 we 
define two pairs of canonical variables for the Ishimori model, re- write the model and the 
Hamiltonian in these variables, and calculate the Hamiltonian on some of the simplest 
known solutions. This is preceded by a discussion of the physical interpretation of the 
model. The Appendix contains a Lax pair for an "extended" system of the deformed 
Heisenberg magnet model. 



2. DEFORMED HEISENBERG MAGNET EQUATION 

a). Canonical variables. 

Let us consider the deformed Heisenberg magnet equation [7|l]: 

Si = SAS,.,. + -SAS,. (2.1) 

X 



Here x = \/x\ + x\ > 0, xi, X2 are the Cartesian coordinates on the plane, S{x,t) = 
{Si, S2, S3), |S| = 1. 

The phase space for this equation is generated by initial data (5*1, 5*2, S'3) subject to 
the constraint |S| = 1. The Poisson brackets of the canonical variables Si in the model 
satisfy the standard relations: 

{Si{x),Sj{y)} = -eijkSk{x)5{x - y), i,j,k = 1,2.3, (2.2) 

where eijk is the fully antisymmetric third rank tensor. For any two functionals F, G we 
then have 



^This equation can be thought of as a cyHndrical-symmetric reduction of the (2+l)-dimensional non- 
integrable Landau - Lifshits equation, St = S A AS. The relation between the latter and the system of 
coupled nonlinear Schrodinger equations in the dimension (2+1) has been discussed in detail in [17]. 
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On taking into account (2.2)-(2.3), one can represent the equation (2.1) in the following 
Hamiltonian form: 

St = -{H,S}, (2.4) 

where the Hamiltonian H is given by 

1 1"°° 

H = - xSl dx. (2.5) 

2 Jo 

Let us now define a new dependent complex-valued variable, 

w{x,t) = , (2-6) 

which is, at each fixed moment of time t, the stereographical projection of the unit sphere 
onto the complex plane, w : §^ ^ C U {oo}. 

In terms of this variable the equation (2.1) can be rewritten as 

w'^w 1 

iwt = w^^ - 2 + -Wx, (2.7) 
1 + \w\^ X 

and the Poisson brackets corresponding to (2.2) take the fornix 

{w{x),w{y)} = {w{x),w{y)} = 0, {w{x),w{y)} = --{1 + \w\^f6{x - y). (2.8) 
The bracket (2.3) then becomes 

{F,G} = --j Ml + l-(.-)l V - ^^^] . (2.9) 
and the evolution of the system will be described by the equation 



with the Hamiltonian 

H = 2 X ] ] dx. 2.11 

Jo (l + |w;|^)^ 

It should be noticed that the following "complex extension" of the system (2.1) is of 
interest of its own 

Iris 1 . _ Islr 1 

iTi — r^x — - — ; 1 — i^x, ist — —Sxx + :; — ; Sx- (2-12) 

1 + rs X 1 + rs X 

From (2.8) we obtain the Poisson brackets of variables r h s in the form 



^ In the derivation of (2.8) we use the relations {S±{x), S±{y)} — 0, {S+{x), Sz{y)} — —iS+{x)S{x — 
y), {S+{x), S-{y)} — 2iS^{x)5(x — y), where S'± = 5*1 ± iS2, and the Leibnits's rule. 

^In absence of the nonlinear component the second equation in (2.12) can be interpreted as the free 
Shrodinger equation with an effective mass. It is evident then that the first equation can be obtained 
from the second by complex conjugation. 
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{rix), siy)} = -1(1 + rsfSix - y), {fix), siy)} = i(l + rsfSix - y). (2.13) 



The system (2.12), as well as equation (2.7), is completely integrable (see Appendix) and 
have a Hamiltonian structure with the Hamiltonian 



X- 



1 + rs)'- 



dx 



and the equations of motion 



n = -{H, r}, St = --{H, s}, 

X X 



(2.14) 



(2.15) 



and can be considered a model of the system of two coupled deformed Heisenberg's mag- 
nets. 

The Poisson brackets (2.13) can be found from the expression for symplectic two-form, 



$ = i 



dr A ds 



dr A ds 



;i + rs)2 (l + fs)2) 



dx, $ = d(f, 



where 



-'f 

Jo 



ds 



ds 



dx, 



(2.16) 



(2.17) 



-s(l-l-rs) s(l-|-rs))J 

thus (2.16) and (2.17) agree with the corresponding expressions obtained in [18] for the 
standard Heisenberg magnet. 

Notice also that the equation (2.7) is a bi-Hamiltonian system: 



SHi 
5w 

SHi 
SiS 



G 



SH2 
■ , 
2 \ SH2 I ' 



(2.18) 



where Hi coincides with H given by (2.11), the second Hamiltonian H2 reads as 



H2^-i / X dx, (2.19 

and Gi = G\{w,w), G2 = G2{w,w) are the so-called Hamiltonian operators of the form 



-1 



(2.20) 



x{i + \w\^y\-i 

An expression for the matrix operator G2 can be obtained from results in paper [19] 
on the standard Heisenberg magnet but is too cumbersome to be written here. Let us 
just mention that its matrix entries contain a differential and an integral operator thus 
rendering it non-local. 

The relations (2.18)-(2.20) mean that the recursion operator of the equation (2.7) 
under the assumption that det G2 7^ is represented in the form 



R — GiGo ^ ■ 



(2.21) 
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Since (2.7) is a completely integrable system, it admits infinitely many integrals of motion, 
{In}'^=i [7], in involution, that is, satisfying {Ij,Ik} = 0. In turn, this allows to obtain 
hierarchies of the Poisson structures, 

/„ = RIn-i, (2.22) 
and the higher equations of the deformed Heisenberg magnet (j = 0, 1, to = ^), 

twt^ = WG2^-^. (2.23) 

ow 

b). Stability of certain solutions of equation (2.7). 

The problem of stability of stationary solutions of the equation (2.7) is of interest since 
the equation contains the independent variable x explicitly. To analyze it, let w = Wgt + w. 
On linearizing (2.7), first on the trivial background Wst = 0, which corresponds, in terms 
of the magnetization vector, to the vector S = (0, 0, 1), we obtain: 

1 . 

iwt{x, t) = w^^{x, t) + -Wx{x, t). (2.24) 

Suppose that (x > 0) 

w{x, 0) = woix), w{0, t) = u'i(t). (2.25) 

Then the equation (2.24) can be solved by the Laplace transformation in the t variable 
under the additional assumption that |w(x, t)| < Me^'-^^ with an M > and sq > 0. 
Solving the arising equation and performing the inverse transformation we find: 

pa+ioo 

wix,t) = — eP'[Coix,p))Joiy^x)] dp, (2.26) 
where Jo(.) is the Bessel function, 

C,{x,p) = -z r e QiO ^€ [ r Qi^) dy ] dx, (2.27) 

Jo Jo Mv-^py) 

Q{x) = -2a/^ (In Jo( V'-^a;))^. + ^, 

Re a > 0, the path of integration is any straight line Rep = a > sq > 0, and the integral 
in (2.26) is understood in the sense of the principal value. It is not difficult to see that the 
logarithmic divergencies arising in the exponentials when integrating at the lower limit in 
(2.28), cancel each other. 

It follows from (2.26) that for a fixed x the function \w{x, t) \ grows with the t increase, 
and, as in [7], we obtain that the solution is unstabl^: an arbitrary localized initial 
perturbation of the system can grow indefinitely as the time passes. 

We now proceed to analyze stability of the stationary state Wgt = ie^^^^^ where 9{x) = 
ln(x) + 6*0, ^0 ^ is a constant. This solution is an example of a spiral-logarithmic 



"'Of course, the stability of that Hnearized "non- autonomous" equation is meant. 
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structure found in [10]: S = (sin 6', cos^, 0)0. On linearizing the equation (2.26) on this 
background, we have, 

_ _ 1 Q-i^i^) 

iwt{x,t) = Wxx{x,t) + -Wx - i 5—. (2.29) 

X x^ 

This equation only differs from (2.24) by the presence of a non-homogeneous term. Hence, 
its general solution is a sum of (2.26) and a partial solution. It follows that it will be 
unstable as well. 

Notice then, that the equation (2.7) admits a solution periodic in t of the form w{x, t) = 
W{x)e^''^ with k a real constant, provided that the equatior@ 

2W'^W 1 

- 1 , rwi2 + -W. + kW = (2.30) 
1 + |VV p X 

has a solution. This suggests that the study of the linearized stability is insufficient. 
The analysis of the nonlinear stability requires more subtle methods [see, e. g. [20] and 
literature cited therein]. 



3. LANDAU-LIFSHITS MAGNET 

a). Canonical variables. 

The fully anisotropic model of Landau-Lifshits has the form 

St = S A S^^ + S A JS, (3.1) 

where J = diag{Ji, J2, J3) are diagonal 3x3 matrices, and Ji, J2, J3 are parameters of 
the anisotropy, Ji < J2 < J3. 

The Hamiltonian for (3.1) can be written in the form, 

1 

H=- {Si- SJS)dx, (3.2) 



2_ 

or, using the variable w defined in (2.6), as0 

H= r ^"'"'.tl' - (3.3) 

where 



2(kl 




+ w^) — 7 w 




(1 + 


\w\ 


2^ 



^ Using (2.11), it is easy to check that the Hamiltonian logarithmically diverges on this solution in 
both limits and, thus, requires a regularization. 

^Removing the nonlinear term we obtain here the stationary Schrodinger equation with the Coulomb 
potential and an effective mass. 

''It is well-known [21], that this model is one of the most general completely integrable models admitting 
2x2 -matrix Lax representations. 

^We assume here that w is a slowly decreasing function. In the case of a decreasing w one should add 
J3 — 4/3 to the density of the Hamiltonian. 
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J2 — Jln J3 J J1 + J2 /g ^\ 

« = — ^ — ,P = ^,1 = J3 ^ — • (3.4) 

Taking into account (2.9), from this we obtain the following equation of motion for 
Landau-Lifshits magnet model 

zwt = t{H, w} = -^(1 + kH'^, (3.5) 
2 ow 

ovE 

iwt = w^^ - 2 , , I |„ 'yw. (3.6) 

1 + 

Let us consider an implication of this form of the equation. Obvious transformations 
lead to the following relation which contains the parameter a only, 



t{\w\')t = {w,w - WW,), + 2 ^'f , ."^'^'^ + 2aiw' - w'). (3.7) 

1 + \wr 



Letting w = pe**^, r;i;e p = p{x,t), ip = p{x,t), p, (p eM, we obtain: 



{p% = 2{p'4>.). - + W sin 2^. (3. 



^P^Px4> x 

Defining the variables R = p^ n Q = 2p'^(j),, we now find the following "conservation 
law": 



Rt = Qx- 2Q[ln(l + R)], + AaR sin( / %) dx. (3.9) 

J -00 R 

It is especially simple when a = 0, which corresponds to the anisotropy of "the easy plan" 
type. 

In a way similar to (2.18), one can produce a bi-Hamiltonian structure for (3.6) with 
Hi equal to H defined by (3.3), the Hamiltonian 



°° WtW — WrW 



^2= / nil 12^1 12 ^^' (3-10) 
, (1 + |iyp)|ti?|^ 



and the Hamiltonian operators 

(1 + HT (-1 0) ^^'^^^ 

and G2 being a matrix integro-differential operator [22]. In terms of the variables w m w 
the recursion operator can be written as follows: 

Y-l 



R = GiG^\ (3.12) 



^Notice that, as well as in the case of the deformed Heisenberg magnet, we are able to obtain the 
corresponding complex extension (see, [18]); we are not going to dwell on that here. 

^°0n taking the complex conjugated equation and neglecting the nonhnear component, one can obtain 
the nonstationary Shrodinger equation with the potential V = —7 = const. 
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It produces an hierarchy of the Poisson structures similar to (2.22) and higher Landau- 
Lifshits equations similar to (3.5). 

b). The Dispersion relation. Stationary Landau- Lif shits equation. 

Linearizing the equation complex conjugate to (3.6) and choosing w = w{x,t) in the 
form w ~ exp{i{kx — ut)}, we have, 

u = k^--f, (3.13) 

which gives a dispersion relation for the Landau-Lifshits equation which is typical for 
magnets with an exchange interaction [23]. In our case the group and phase velocities are 
given by Vg = du/dk = 2k, Vph = oo/k = k — ^y/k, respectively (the latter is infinite for k = 
0), implying that there is a dispersion in the system. The propagation of a magnetization 
wave in this model is possible under the condition fc^ > 7 = J3 — ( Ji + J2)/2 > 0. 

Letting w = w{x — fit) = w{^) in (3.6) , where fi = const is the velocity of a stationary 
profile wave, we obtain the equatioiuj 

w{w] + a) — aw^ — -yw 

wcc + i^Wf — 2 ^ j — jw = 0. (3-14) 

1 + \wf 

From this it is not difficult to obtain that 

2 —2 —2 2 

{w^w - w^w)^ + ifx (l^n^ - 2 ^^|^|2 — + 2a(?^ - w^) = 0. (3.15) 
Letting w{^) = p e*"^, where p = p{^), (p = 0(0^ P ^ G IR, we then have: 

2(p2),0, + 2p20 + ^(p2) _ g^Pi^ + 4a sin 20 = 0. (3.16) 

1 + p"^ 

Let /i 7^ 0. Then, obviously, p = const, = 7m/2, n = 0, ±1, ±2, ±3, satisfy (3.16). For 
= 00 = const we obtain, p^ = Pq — (4a//i) sin(20o)^, where po = const; npn p = po = 
const (3.16) is reduced to the equation of the pendulum: 0^^ + {2a/ po) sin(20) = (the 
existence of other solutions remains an open problem). 
Let now p = 0, then from (3.16) it follows that 

^-r^ = Ci, 0^c + Ci05 + 2asin20 = O, (3.17) 
p^ 1 + p^ 

where Ci is arbitrary constant. The first of these equations can easily be integrated: 

pi 2(^) = 1(1 ± ^1 _ 4e-2(Ci5+C2) ) e^i«+^2, (3.I8) 

where C2 is another arbitrary constant (we assume that e"^^*-^!^"*"'-^^) < 1/4), and the second 
equation, which coincides with the one of the pendulum with the friction 0, admits, in 
particular, solutions of the form = 7m/2, n = 0, ±1, ±2, ±3. 

Thus, solutions of the stationary Landau - Lifshits equation have fairly non-trivial 
structure in the generic (fully anisotropic) case. Their further study could bring a solution 



Stationary equations of another form for the Landau-Lifshits hierarchy were considered from the 
viewpoint of the Lie-algebraic approach in [24] . 

^^In the partial case Ci = this equation, can obviously be integrated in terms of the eUiptic functions. 
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to an important problem in the theory of dynamical systems - that of construction of the 
phase graph for the equations (3.14) and (3.6)0. The same applies to the deformed 
Heisenberg magnet from the previous section. 



4. ISHIMORI MAGNET 

a). Physical and geometrical interpretations. 

The Ishimori magnet model in terms of the magnetization vector has the form: 

= S A {S^x + Ct^^yy) + UySx + UxSy, (4.1) 

- a^Mj/s/ = -2a^S(Sx. A Sj/), (4.2) 

where S{x,y,t) = {81,82,33) is a three dimensional vector, |S| = 1, u = u{x,y,t) is an 
auxiliary scalar real- valued field, and the parameter takes values ±1. The system is 
called the Ishimori-I magnet (MI-I) in the case = 1, the Ishimori-II magnet (MI-II) in 
the case = —1. Mathematically, each of these cases corresponds to different types of 
the equations (4.1) and (4.2). 

The topological charge of the model (4.1)-(4.2), 

Qt = ^ [ [ S{S., A Sy)dxdy, (4.3) 
4vr Jr2 J 

is invariant under the evolution of the system. Since the homotopy group of the unit 
2-sphere 712(6'^) coincides with the group Z of integers, the number Qt must be integer. 
According to (4.3), the scalar function u = u{x,y,t) is related to the density of the topo- 
logical charge production. The derivatives u^, Uy in (4.1) play role of friction coefficients. 
Thus, (4.1) can be interpreted as an equation of forced (by the friction power) precession 
of the magnetization vector, and the system (4.1)-(4.2) is self-consistent. 

From the physical viewpoint, it is easy to see that there is a non-local interaction in 
this system, on top of a local (exchange) one. The mechanism of the former is unclear. 
Nevertheless, the study of such systems is justified since stable localized two-dimensional 
magnetic structures are observed in experiments. An argument in favor of this assertion 
is the above-mentioned gauge equivalence of the MI-II model and the DS-II model, which 
describes quasi-monochromatic waves on the fiuid surface [5], and also a link found in [26] 
between the MI-I model and the nonlinear Schrodinger equation with magnetic field. 

Also helpful is another, hydrodynamical, interpretation of the model (4.1)-(4.2). Namely, 
let Uy = —Vi, Ux = V2, hence w{x,y) = {vi, V2) is the velocity field of a fiuid. Then the 
MI model can be rewritten as follows: 

St + ViSx - V2Sy = S A {Sxx + oi^^yy), 

(4.4) 

^^Phase graphs of the equation (3.1) in the case of partial anisotropy have been studied in [25]. Phase 
graphs in the fully anisotropic case have apparently not been considered yet. 
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V2x + Oi^viy = -2a^S{Sx A Sy). 

If we define the stream function of the flow, vi = — Xiy, V2 — Xix, then the equation (4.2) 
with — —1 (the MI-II model) implies the Poisson equation 

Xixx + Xij/2/ = 2S(Sa; A Sy), (4.5) 

that is, the stationary (the time t is a parameter here) vorticity equation with a source 
in the right hand side of the magnitude proportional to the density of the topological 
charge production (details on the equation of planar hydrodynamical vortex can be found 
in [27]).^ 

Let F{x, y, t) = 28(8^: A Sj^). On taking one of the expressions of the form ±6=*=^^, e^^ — 
g-2xi^ isinhxi, icoshxi, isinxi, icosxi, tov F{x,y,t), we obtain a closed completely 
integrable equation of elliptic type for the function Xi- The solution of an appropriate 
boundary-value problem for this equation must satisfy the additional condition 

-/ [ Axi{x,y)dxdy^No, No e Z, (4.6a) 



or 



(r = ^/ x"^ + y"^) 



lim ^ / {Xixdy - Xiydx) = Nq. (4.66) 
b). New canonical variables. 

Let us now consider another canonical variables. First, we pass from the variable 8 to 
new variables p h g (p, g G M) in (4.1)-(4.2), setting [28]: 



Ss{x,y,t)^p{x,y,t), S+{x,y,t) ^ ,Jl - p\x,y,t) e^''^^'^'*). (4.7) 

Expressions for Poisson brackets of the quantities p n q follow directly from (2.2), on 
taking into account that the problem is two-dimensional, 

{p(r),g(r')} = (5(r-r'), {p(r),p(r')} = {g(r), g(r')} = 0, v^{x,y), (4.8) 
and then for the any two functionals F and G one can obtain: 

rr. ^1 f f, 5F 5G 5F 5G , , , , 
J^2j 5p{r)5q{r) 5q{r) Sp{r) 

In terms of this variables the MI model (4.1)-(4.2) can be rewritten as a Hamiltonian 
system. 



6H Pxx + a^Pyy P {pi + o'pI) . 2 ^ 2 2N ^ 

g* = ^ = f3^2 (i_p2)2 - P + « %) + ^y<lx + u^qy. 

5H 

Pt = -— = (1 - p'^){qxx + oi^qyy) - 2p {Pxqx + a^Py^y) + UyPx + UxPy, (4.10) 

^xx ^ ^yy 2*3; {Pyqx Pxqy)i 
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and for the topological charge we will have: 



Qt = ^ [ [ {Pyqx - PxQy) dx dy . (4.11) 
4vr Jk2 J 

Here the Hamiltonian H has the form : 

H = H, + H,, H, = U [ fl±^ + (1 _ p^)^ql + a\l)]dxdy, 

(4.12) 

H2 = - I I [a^A^ + B^] dx dy, 

where A = Ux, B = —a^Uy, so that A^ + By = 2a'^{j)x<iy ~ PyQx)', in this case it can be 
take place the conditions: 



^ = C5y{y - y')6{x - x), ^ = D5{y - y')5x{x - x'), 



^ = E5y{y - y')6{x - x'), ^ = F5{y - y')5x{x - x'), 



(4.13) 



where C, D, E, F are some functions. Letting D = C, F = E and taking into account 
(4.10), we obtain the following relations on the functions C = C{x, y, t) and E = E{x, y, t) 
(the symbol <, > refers to the scalar product in R^, and T stands for the transposition): 



1 „ 1 

from this we find: 



(4.14) 



C(a;, t) = Cq{u{x, 2/)) + ^ / (^x^jy + Uyqx)ds, 



(4.15) 



^(a;, 2/, t) = Eo{u{x, y)) + Tr^ / (mxPj, + UyPx)ds 



2a2 

where Co, Eq are arbitrary functionals, and the integration goes along the characteristic 
s of the equations (4.14). Assuming that Eq = Cq, we see that the functional Cq must 
obey an additional condition: 



^''Paper [29] contains an expression for the Hamiltonian of the so-called modified MI different from 
(4.1) by the sign in the last but one term. Thus, the Hamiltonian for the model (4.1)-(4.2) seems to have 
been obtained here for the first time, both in the q, p and w, w variables, the latter being defined below. 
Also, in contrast with the modified model, it is easy to see that it is impossible to define the Clebsch 
variables in our case. 
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6Co ,6u 6u 2 

Since u^y 7^ in the generic case, from this it follows that one more condition is necessary: 
6Co/6u 7^ (if, of course, at this 6u/6q ^ 6u/6p). 

Now let us pass to the variable w in (4.1)-(4.2), defined in (2.6) (assuming that w = 
w{x,y,t))^ : 



IWt = Wccx +a Wyy - 2 — — pr-^ + t[UxWy + UyW^), 



(4.17) 



_ 2 _ 2 '^xWy - WxWy 

(1 + \wr)^ 



Then for the topological charge we obtain 

i f fWxWy-WxWy 

Qt = -7- —-^—-^dxdy. 4.18 

2vriiK2i (l + |wP)2 

The non-vanishing of the Poisson bracket for the canonical variables w{x,y),w{x,y) 
comes along as in (2.8): 

{w{x, y),w{x', y')} = -^(1 + HY^i^ - r'), r = (x, y). (4.19) 
This allows to rewrite (4.1)-(4.2) in a transparently Hamiltonian form: 

zw, = -hl + \w\r^-^. (4.20) 
2 ow 



Here H is the Hamiltonian of the form 



WxWx + a WyWy 



H = H, + H,, H, = 2 I I " dxdy, 



H2 = ^ 1 f \a^u\ + mJ] dx dy, 



(4.21) 



and we assume in the course of the derivation of the equations for the model that the 
following conditions, analogous to (4.13), are satisfied: 



Clearly, all three representations of the MI model, (4.1)-(4.2), (4.10) n (4.17) are 
equivalent. 



^^The reflection [w, w) — > (p, q) can by given by relations q — — arctan(i(u' — w)/{w + w)), p 

(\w\^-i)/ii+\wn 
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Notice also that, one can define a "complex extension" of the system (4.17) analogous 
to the ones above. Letting formally wi = w, one obtains, 



2 + 

IWt = Wxcc + a Wyy- 2 -r^ + + M^Wx;, 

(1 + tytyi)^ V y y / 



2 "^("^L + '^^'"^Ih) 

(wii:^ + a wiyy) + 2 — _^ ^2 h i{u^wiy + (4.23) 



_ 2 _/|- 2 ^x'^iy ~ WlxWy 
Uxx " Myj/ - 4m ^ u;tt)^)2 



This system can be interpreted as a model of two coupled Ishimori magnets. Nontrivial 
Poisson brackets follow from (4.19): 



{w{x,y),wi{x\y')} = --{l+wwif6{r-r'), {w{x,y),Wi{x',y')} = -{l + wwif6{r-r'), 

(4.24) 



and the "topological charge" of this model is0 



Qt = ~ / "-"'■:'"'\f' ^:.^y. (4.25) 



The equations of motion (4.23) are Hamiltonian: 



1 , . 1 , \2^H , , 

iwt = --[I + wwi) - — , iwu = -[l + wwi) (4.26) 

2 owi 2 ow 



where 



H = Hi + H2, Hi = 2 / — ^^dxdy, 

H2 = ^ j j [a^ul + M^] dx dy, 



and we suppose that 



6Ux AiWy 6Uy 4:iwx 



Swi a^(l+wtyi)^' 6wi (1 + wwi)^' 

(Jtii; AiWly 6Uy AilWlx 



(4.27) 



(4.27a) 



6w a^(l + wwi)^' 5w (1 + wwi)^ 

Returning to (4.17), we introduce the complex coordinates z = x + iy, z = x — iy, so 
that dz = l/2{dx — idy), ds = l/2{dx + idy), dx dy = {i/2)dz A dz and rewrite (4.18) and 
(4.21) in terms of these variables. 

i). Let = 1, that is, the MI-I model is considered. In this case we obtain: 



^^In general, wi ^ w, thus Qt can be non-integer and even complex. Such a situation, including an 
interpretation of the quantity Qt, should be considered separately. 
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WzWz 

IWt = AWzz - 8 W - 2{WzUz - WzUz), 

1 + \wr 



_ WzWz - WgWz 
Uzz ~r Uzz ^" 



;i + \w 

The topological charge is given by 



2\2 



and the Hamiltonian is 




iWt = 2{Wzz + Wzz) - 4 ^ I W - 2{WzUz - WzUz), 

1 + \wr 



_ ^ WzWz - WzWz 



(1 + kP) 



(4.28) 



^ ill WzWz — Wf'U^z , , . N 



H = 2t[ /' ^^^^- + ^^^^ dzAdz+- I luzU-zdzAdz. (4.30) 

J J (1 + \w\^y 2 

ii). Let = -1. We then have the MI-II model, 

,wt + w'^ 



(4.31) 



The expression for the topological charge coincides with (4.29), and for the Hamiltonian 
we have: 

H = 2tJ J dzAdz-'-J j{ul + n^) dz A dz. (4.32) 

Notice also that the Hamiltonian of the MI-I magnet and its topological charge are 
related, as follows from (4.29) and (4.30), by the inequality of Bogomol'nyi, which is a 
lower estimate for the Hamiltonian taking into account all the dynamical configurations. 
Namely, 



H > AtiQt. (4.33) 

Comparing (4.29) and (4.32) one can see that for the MI-II model such an estimate does 
not exist. 

c). Hamiltonians and topological charges for some of the simplest solutions. 
Equations (4.1)-(4.2) can be interpreted as the compatibility conditions for the follow- 
ing overdetermined matrix systems on the function \1/ = \l/(x, ?/,t): 

= -S^^, (4.34) 
a 

= -2iS^,, + Q^,, (4.35) 

Similar calculations were given in [29] for the case of the modified MI and certain other systems. 
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where Q = Uyl + a^UxS + iaSyS — iS^, = '^{x,y,t) G Mat(2,C), S = ^^=1 5'j(Tj, cTj are 
the standard Pauli matrices, / is the unit 2x2 matrix. By the definition, the S matrices 
have the properties: S = S*, S'^ = I, det S = —1, Sp S = (the asterisk stands for the 
Hermitian conjugation). 

For a future reference, let us provide an expression for the S matrices in terms of the 
w variable: 



( |w| -1 2iD \ 

l+|wP 1+|«)|2 / 

and consider some of the simplest examples of calculations. 

1. Let = —1 in (4.17), that is, the MI-II model is considered. As was shown in 
[26], the conditions 

w, = 0, w, = 0, (4.37) 

are then compatible with (4.31). The first and second conditions here mean the presence 
of instanton and anti-instanton sectors, respectively, in the MI-II model. Consider the 
instanton sector assuming that w{z) = {{z — zo)/A)" [300, where n E Z+, A G C (the Zq 
and A characterize, respectively, the position and size of the instanton). A calculation by 
the relation (4.39) gives: 

Qt = ^/ / (^[""i^pp dzAd-z = n, (4.38) 

and i^i = by (4.32). To find H2, one has to know the function u. Using the second 
equation in (4.31) with Wz = and returning to the Cartesian coordinates, we obtain 

{zo = xo + iyo), 



[(x - xqY + jy - yof^^ 
[A2- + [(x - xo)^ + (1/ - yo) W 



A« = -8y \.J) i;/^..' :.\2i..2 > (4-39) 

which implies that [^^' 



where Gq{x, y) = (l/27r) ln(a;^ + y"^) is the Green function of the two-dimensional Laplace 
operator. Thus, the energy of the instanton solution on the formal level is given by 




(4.41) 



The whole instanton sector is then split into disjoint classes each corresponding to the 
relevant value of the Qt quantity. 

^^The choice of a more general solution, say, in the form of the Belavin-Polyakov instanton (linear- 
fractional function with complex poles) [2], unfortunately, significantly complicates the calculations. 

^^The following integral can be simplified by a change of variables and subsequent contour integration, 
but the remaining integral, apparently, cannot be calculated explicitly. 

Obviously, the Hamiltonian is positive in the domain where \uy \ > \ux\. 
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2. Let us consider the MI-I model (a^ = 1) and show that instanton solutions exist in 
there as weli£l|@. Indeed, for Wz = the system (4.28) is reduced to the following one, 

wt = -2uzWz, Uzz + Uzz = 4 / ^ . (4.42) 

(1 + 

Differentiating the first relation in z, we obtain that Uzz = 0, whence the compatibility is 
achieved if 

= \, !f m2 > (4-43) 

(1 + \w\'^y 

or 



u{x,y,t)=A dx' dx" + foiy,t)x + c,, (AAA) 





\w 




2 


(1 + 


\w\ 





where Ci is an arbitrary constant, and /o(., .) is an arbitrary function. For the instanton 
solution w{z) of the same form as in the previous case the number Qt = and the 
function 



An' r , , r' „ \{x" - x,f ^ {y - y^ff-' , , 
^ J-oo J-oo + - xo\' + \y - yol'ry 

The expression for the Hamiltonian takes the form, 

H = A7in + \ I l{ul + ul)dxdy. (4.46) 
A Jr2J 

3. Let us calculate the topological charge and the Hamiltonian for a solution of the 
form of a spiral structure. Namely, let S = (0, sin$i, cos$i), where $i = 6ot + aox + 
Pay + 70, C(o, Po, lo, So & M. are parameters, that is, the solution is a two-dimensional 
spiral structure [16]; then, according to (2.6) (see also (4.36)), w{z,z) = itan($i/2), $1 = 
dot + az + az + 'Jo, a = ao/2 + (3o/{2i). It follows from (4.29) that Qt = 0. 

To determine the function u = u{x, y, t), one has to substitute the function w{z, z)) 
in the equations (4.28), (4.31), which gives two linear equations for u. Assuming their 
compatibility and integrating, we find (a^ = — 1)[16]: 

uix, y) = goiy + —x) + f gi{y{s') + —x{s'),t) ds\ (4.47) 
ao Js ao 

where go, gi are arbitrary functions such that go is constant on the characteristic y + 
{l3o/oLo)x = const, and s is the characteristic taken to be the integration path. 
Similarly, for = 1 we have: 

^^This is not surprising, albeit apparently went unnoticed in the literature, given that in the "static 
limit" the MI-I model turns into the elliptic version of the nonlinear 0(3) tr- model for which the instanton 
solutions were constructed initially. Notice also that the model was solved by the inverse scattering 
method in [31]-[32]. 

^^Prom the viewpoint of the higher-dimensional inverse scattering method and the dressing procedures 
for solutions, the characteristic variables ^ = {y — x)/2 and ij = {x + y)/2 [14], [15] are more natural than 
z and z. 
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u{x, y) = g,{y - ^x) + [ gsiy{s') - t) ds', (4.48) 

where g2, gs are arbitrary functions, and g2 is constant on the characteristic si y — 
{Po/ao)x = const. 

The substitution w = w{z, z) in (4.30) and (4.32) (in both MI-I and MI-II cases) 
leads to divergence of the Hamiltonian Hi, and, therefore, that of the Hamiltonian H as 
a whole, since the functional H2 is finite. 

4. As was first shown in [33] (see also [14], [16]), in the reflectionless section of the 
MI-II model the system (4.34)-(4.35) can be written in the form 

^, = 0, ^i + 2z^,, = 0. (4.49) 

In turn, the latter system has well-known polynomial solutions describing vortex states, 
(^ = = 1,2, ^22 = ^11, ^^12 = -^2i)[33]: 

j=0 m+2n=j 

(4.50) 

Ml /) 1 1 

^-^ ^-^ m\n\ 2 2 

jr=0 m+2n=j 

where Ni is an integer. Mi = Ni — 1, aj, bj are complex numbers, and the inner summations 
run over all m, n > such that m + 2n = j. In particular, in this simplest case Ni = 1 
it follows that 



^'n = ao + aiz, ^21 = &o- (4.51) 

Let us employ now a dressing (say, the Darboux dressing [16]) relation for the matrix 
8 0,5"= ^S"^^^^"^, where S^^^ is the initial solution of the system (4.1)-(4.2), assuming 
that S^^^ = (T3. This leads to the so-called "one-lump" stationary solution, which we 
write here in terms of the stereographic projection, 

(4.52) 



ai{z - Zq) + do 

where do — CLq — fli^o, Zq is the coordinate of the vortex center on the complex plan. 

It is known [33], that Qt = 1 for such a solution. Let us calculate the function 
u = u{x,y). In the case under consideration the second equation in (4.31) is reduced to 
the following, 



2| 


\ai\ 


2 


1^0 


2 


{\ai{z - Zq) + do\ 


2 + 


l^ol 


2^2 



= tt-t:. — z^—TTT-in^- (4-53) 

From this we find 



^^Its structure is identical, at least, for all the models of magnets treated here and all known methods 
of their solutions. 
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«(a;,|/) = 2|ain6or / / G^{x - x' ,y - y')j^^, ^L^^^^ dx' dy' . (4.54) 

Taking into account that Hi — 0, we now obtain from (4.32) the Hamiltonian of the 
one-lump solution in the form 

f [{ul- ul) dx dy. (4.55) 

The left hand side is positive in the planar domain where \uy\ > \ux\ in the same way as 
in (4.41). 



5. CONCLUSION 

The results for the Ishimori magnet show, in particular, that the Hamiltonians and 
topological charges cannot always be calculated analytically even for the simplest so- 
lutions, and numerics are required for specific Cauchy problems. In this respect, it is 
especially interesting, in our opinion, to check the hypothesis [16] of possibility of a phase 
transition in the model which involves a change of topology and symmetry properties of 
the system. 

Concerning the "extended" systems (2.12) and (4.23), we would like to point out that 
if the initial systems are gauge equivalent to the nonlinear Schrodinger equation with an 
integral nonlinearity [8,9] (the Davy-Stewartson-II model [5] in the IM-II case), then it is 
interesting and important to find objects gauge equivalent to the extended systems. 

Overall, the representations considered in this paper can, hopefully, be useful in study- 
ing other (1+1)— (and more realistic (2+1)— and (3+1)—) dimensional models of mag- 
nets, (T-models and chiral fields, including nonintegrable cases. 

The author is indebted to P. Kulish for support. 

APPENDIX 

We provide here an expression for a Lax pair for the system (2.12). First, define the 
function D = D{x, t) e Mat{2, C), x > 0, of the form 

/ sr-l 2s \ 

D(x,t) = (]tr ■ (Ai) 

V l+sr 1+sr / 

The matrix D has the following properties: = I, Sp D = 0, det D = —1. However, 
unlike the matrix S, it is not Hermitian. A straightforward calculation shows that (2.12) 
is the compatibility condition for the following overdetermined linear system of equations: 
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*, = -^AD*, *,= (^D,D + |a2D)*, (A2) 

where "if — ^(x, t, A) G Mai{2, C), A = A(x, e C. This means that the parameter 

u plays the role of a "hidden" spectral parameter, so that the conditions 

A ^ 2A^ , . 

A. = -, At = , ^.3 

X X 

or 



2(^ + 1^)' 

are fulfilled, and the matrices D satisfy the equation 

iDt = l-[D,D,,]--D,D. (A5) 
z X 

Thus, we have a non-isospectral deformation of the associated linear system (the case 
of a single deformed Heisenberg magnet was considered in [8]). Notice also that so far 
the inverse scattering transform method has not been applied to such systems at the full 
scale. 



(A.4) 
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